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DAMPING  ANALYSIS 
SINGLE  AND 


OF  WAVES  IN  FLUID-LOADED 
CONSTRAINED  PLATES 


1.  INTRODUCTION 

The  technique  of  constrained-layer  damping  is  based  on  the  addition  of  two 
layers  to  a  base  plate  or  bar:  the  first  layer  consists  of  a  viscoelastic 
material  with  high  loss  factor  in  its  shear  modulus;  the  second  layer  is  made 
of  a  stiff  material  that  increases  the  shear  in  the  viscoelastic  layer,  and 
thus  promotes  viscoelastic  dissipation  (Fig.  1). 

An  early  analysis  of  constrained-layer  damping  was  developed  by 
E.  M.  Kerwin  in  1959.  [1]  (A  review  of  this  technique  is  found  in  Ref.  2.)  It 
is  limited  to  thin-plate  theory  and  does  not  include  extensional  waves. 
Recently,  a  model  was  developed  by  this  author [3]  that  combines  a  description 
of  the  base  plate  by  exact  elasticity  theory  with  the  Kerwin  treatment  of  the 
damping  mechanism.  The  hybrid  model  lends  itself  quite  readily  to  inclusion 
of  the  effect  of  fluid  loading.  It  is  often  tacitly  assumed  that  the 
radiation  and  viscoelastic  damping  in  a  fluid-loaded  constrained  plate  are  the 
same  as,  respectively,  the  radiation  by  a  single  plate  in  the  fluid  and  the 
viscoelastic  damping  of  a  composite  plate  in  vacuum.  This  ignores  the 
interaction  of  fluid  and  composite  plate,  with  possible  consequences  for  the 
computation  of  attenuation  and  damping  efficiency. 

The  present  study  develops  a  model  for  a  hybrid  constrained  layer  that  is 
coupled  to  a  fluid  loading  the  plate  on  either  or  both  faces.  The  model  is 
formulated  and  applied  to  the  attenuation  of  flexural  waves  in  two  cases:  on 
one  hand,  the  combination  of  the  results  for  a  fluid-loaded  single  plate  with 
those  for  a  composite  plate  in  vacuum;  on  the  other  hand,  the  results  for  a 
fluid-loaded  composite  plate. 

While  studying  the  attenuation  of  wave  propagation  in  a  constrained  layer 
under  fluid  loading,  it  was  found  that  at  low  frequency  (well  below  the 
coincidence  frequency)  the  attenuation  coefficient  again  reaches  high  values. 
There  does  not  appear  to  be  a  wide-spread  awareness  of  this  phenomenon  (which 
is  also  present  for  a  single  plate  in  a  fluid) ,  although  it  is  found  in  the 
literature.  Therefore,  discussion  of  the  fluid-loaded  hybrid  model  is 
preceded  by  a  discussion  of  the  radiation  by  a  single  plate  into  a  fluid 
according  to  thin-plate  and  exact-elasticity  theories.  This  may  also  function 
as  an  introduction  to  the  literature  on  the  study  of  fluid-loaded  thin  plates. 

2.  FLEXURAL  WAVES  IN  A  FLUID-LOADED  INFINITE  THIN  PLATE 
2 . 1  Theory 

The  theory  of  the  effect  of  fluid  loading  on  the  propagation  of  flexural 
waves  in  an  infinite  plate  in  the  thin-plate  approximation  is  extensively 
described  in  the  literature.  A  general  discussion  of  the  problem  may  be  found 
in  Refs.  4  and  5.  More  detailed  studies  of  the  transfer  function  describing 
the  forcing  of  fluid-loaded  plates  are  given  in  Refs.  6  through  11. 
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For  a  discussion  of  the  basic  theory,  the  approach  presented  in  Ref  5  is 
followed  here  A  straight-crested  flexural  wave  is  traveling  in  the  x- 
direction  in  an  infinite  plate  with  thickness  h,  density  /?  ,  Young’s  modulus 

E,  Poisson’s  ratio  and  bending  stiffness  D;  D  =  Eh  / 1 2 ( 1  -  f/).  A  fluid 

with  wave  speed  cq  and  density  p is  adjacent  to  the  plate  at  the  side  of 

positive  z;  z  is  the  coordinate  perpendicular  to  the  face  of  the  plate;  the 
surface  of  the  plate  exposed  to  the  fluid  is  at  z  =  h/2.  The  equation  for  the 
displacement  w  of  the  plate  in  the  z  direction  is  [5] 


DV4w  +  o  h 
rs 


92_w 

at2 


-p(z  =  h/2) , 


(1) 


where  p  is  the  pressure  fluctuation  in  the  fluid. 

It  may  be  noted  here  that  the  assumption  of  one-sided  fluid  loading  is 
incompatible  with  the  exclusion  of  extensional  waves,  but  the  ensuing 
discrepancy  is  small  for  the  frequency  range  in  which  the  thin--plate 
approximation  is  acceptable . 

Continuity  of  normal  displacement  at  the  plate  surface  for  z  =  h/2  leads 
to  the  boundary  condition 


= "  1/po  If  at  z  =  h/2- 


(2) 


If  the  wave  in  the  plate  is  harmonic,  with  time  and  x  dependence  expressed 
by  expi (wt-kx)  (where  k  is  the  wavenumber),  the  corresponding  compressional 
wave  in  the  fluid  is  given  by 


P+  =  PQ+  expi|>t  "  Kxx  -  Kz(z  -  h/2)] ,  (3) 

where  the  subscript  +  refers  to  the  positive  side  of  the  plate,  in 
anticipation  of  a  later  introduction  of  fluid  at  the  negative  side,  indicated 

with  subscript  Here  K  and  K  are  the  components  of  the  wave  vector  ft  in 
the  fluid.  x  z 

Matching  the  normal  displacement  at  the  boundary  in  plate  and  fluid  leads 
to  the  condition  that  k  =  K  .  For  reflection  and  transmission  problems  of 

plane  waves  k  and  ft  are  real,  and  one  may  set  K  =  k  sin#  and  K  =  k  cos#, 

X  o  z  o 

where  kQ  is  the  wavenumber  in  the  fluid,  with  kQ  =  w/cq,  and  #  the  angle  of 
incidence  of  the  incoming  wave.  For  a  (decaying)  free  wave  in  a  plate  k,  ft, 
and  #  are  complex. 
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To  obtain  a  dispersion  relation  with  real  coefficients,  one  introduces  a 
quantity  r  related  to  by  r  =  iK^/k^  [5,10].  The  boundary  condition  Eq.  (2) 

then  gives 


2 

p^U)  w  =  -kQTp+  at  z  =  h/2.  (4) 

This  enables  one  to  replace  p+  in  the  wave  equation  [Eq.  (1)]  resulting  in  the 
dispersion  relation 

k4D  -  ^hw2  =  py  /(kQr) .  (5) 

Since  the  wave  in  the  fluid  satisfies  the  Helmholtz  equation,  it  follows 

that 

K2  +  K2  =  k2  -  k2  r2  =  k2  ;  (6) 

X  z  o  o  }  K  J 

2 

and  thus  one  may  substitute  for  k  in  Eq.  (5) .  This  results  in  an  equation 
for  r  in  the  form 

t5  +  2 r3  +  (l-n"2)r  -  etr3  =  0  ,  (7) 

where  fl  is  a  non-dimensional  frequency  (1  =  wfw  .  The  angular  frequency  UJ  is 

c  c 

that  frequency  for  which  the  wave  speed  c^  of  flexural  waves  in  an  unloaded 
thin  plate  is  equal  to  the  sound  speed  in  the  fluid  (coincidence) .  Since 

cb4  =  w2D/(/>gh)  ,  (8) 

one  has 

“c2  =  P'K*/*  ■  (9) 

The  material  parameter  e  is  given  by 

6  =  Poc o/(Vsh)  •  (10) 

It  is  a  measure  for  the  influence  of  fluid  loading;  for  e  =  0,  the  wave  speed 
c  in  the  plate  reduces  to  c^,  the  wave  speed  in  vacuum. 
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Alternatively  e  may  be  written  as 


e  =  -i-  (p  fp  )  (c  /c  ) 
^2  °  rs  p  cr 


(ii) 


and  thus  e  is  independent  of  the  plate  thickness  (as  long  as  the  plate  may  be 
considered  "thin";  i.e.,  kh/2«l)  Here  c  is  the  wave  speed  for  extensional 

2  o  P 

waves  in  a  thin  plate,  c  =  E/[(l-i/  )p  ]. 

P  s 

2.2  Character  of  the  Roots  of  the  Dispersion  Relation 


The  dispersion  relation  [Eq.  (7)]  has  five  roots,  not  all  of  which  have 
physical  significance.  Since  its  coefficients  are  real  and  the  degree  is  odd, 
there  is  at  least  one  real  root.  The  quantity  T  appears  in  the  spatial 
dependence  of  the  pressure  in  the  fluid  through  the  factor  expf-k^Tz)  and  thus 

only  real  roots  with  positive  value  for  T  are  physically  acceptable, 
representing  an  evanescent  wave  in  the  positive  z  direction.  If  there  are 
complex  roots,  they  have  to  appear  in  pairs  of  complex  conjugates;  and  the 
choice  between  the  two  members  of  a  pair  is  dictated  by  the  fact  that  one 
expects  wave  propagation  in  the  positive  z  half-space  when  there  is  fluid  only 
at  that  side  of  the  plate.  As  a  consequence,  the  imaginary  part  of  r  should 
be  positive.  The  constraint  for  the  real  part  of  r  issues  from  the  following 
observation.  If  one  considers  a  wave  in  the  plate  in  the  positive  x 
direction,  the  complex  wavenumber  k  =  k'  -  ia  is  such  that  k'  and  a  are  both 
positive  to  assure  a  decrease  of  the  wave  amplitude  in  the  x  direction.  Since 


k2  - 


,22 
k  r 
o 


=  kQ,  the  product  of  the  real  and  imaginary  parts  of  k  r  is 

equal  to  -akr  ;  thus,  the  real  part  of  r  has  to  be  negative  This  in  turn 

makes  the  imaginary  part  of  K  positive,  and  the  amplitude  of  the  wave 

z 

increases  with  distance  from  the  plate.  This  may  appear  physically 
unrealistic  but  it  is  not:  Fig.  2  illustrates  the  explanation.  Waves  farther 
away  from  the  plate  originate  in  a  location  at  the  surface  of  the  plate  where 
the  amplitude  of  the  wave  had  a  larger  value 


The  real  root  for  T  represents  a  purely  evanescent  wave  even  for  an 
infinite  plate.  The  complex  root  poses  a  problem  in  interpretation:  it  has 
the  appearance  of  a  radiating  root,  but  only  for  an  infinite  plate  can  this 
name  be  maintained.  For  an  infinite  plate,  the  distinction  between  nearfield 
and  farfield  is  lost,  and  the  "radiating"  wave  is  always  under  the  influence 
of  the  plate.  For  a  finite  plate,  a  wave  with  wave  speed  less  than  the  plane 
wave  speed  cq  cannot  exist  in  the  farfield  and  thus  is  not  truly  a  radiating 

wave.  This  was  pointed  out  by  Pierucci  and  Graham  [12] .  In  the  remainder  of 
this  report  the  term  "radiation"  will  be  used  in  connection  with  the  presence 
of  a  valid  complex  root,  but  it  should  be  realized  that  its  behavior  for  a 
finite  plate  needs  further  scrutiny. 

Setting  Kz  =  +  iK£  one  may  conclude  that  the  angle  5  of  the 

"radiated  waves"  in  the  fluid  with  the  normal  is  given  by 


tan  6  =  k1 /K1 
1  z 


(12) 
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and  the  wave  speed  of  this  wave  in  the  fluid  is 

Cf  =  kcc0/[(k')2  *  (xy2]1/J  .  (13) 

Of  course,  these  statements  are  all  a  consequence  of  the  infinite-plate 
assumption. 

This  discussion  may  be  illustrated  by  the  asymptotic  behavior  of  T  (=  iK  ) 
for  f)+0.  On  physical  grounds  one  would  expect  that  for  low  frequency  the 
inertia  of  the  plate  may  be  ignored  by  dropping  the  second  term  of  Eq.  (1) . 

If,  moreover,  the  compressibility  of  the  fluid  is  ignored,  which  is  equivalent 

to  setting  c/c  «1,  it  follows  that  r2  +  1  -  r2  ,  and  the  polynomial  equation 
^  5  — 3 

for  V  reduces  to  r  =  efl  [5] .  The  roots  of  this  equation  are 


r 


e1  /5  /5e2n*ri/5 


n  =  6,  1,  2,  3,  4. 


(14) 


It  follows  from  the  above  discussion  that  only  the  real  root  =  e1/5f]  3 
and  the  complex  root  T  =  61  ^  ft  3 ^  (cos  144°  +  i  sin  144°)  are  acceptable. 

V 

One  also  infers  that  the  non-dimensional  attenuation  for  the  radiating  wave, 
given  by  a/ k' ,  approaches  tan  36°  =  0.7265  for  fl+0. 

One  could  also  express  the  dispersion  relation  in  terms  of  the  wavenumber 

2  9  2  2  2 

k  in  the  plate.  Since  k  r  =  -K  =  k  -  k  ,  one  can  substitute  for  T  in 

o  z  o 

Eq.  (5)  and  square  both  sides  of  the  equation  to  obtain  a  polynomial  equation 

2 

in  terms  of  rj  =  (k/kQ)  . 

J?5  -  >74  -  2»73n"2  +  2>72ir2  +  J7(T4  -  (e2  +  fi2)fi"6  =  0  (15) 


One  finds  five  roots,  of  which  only  two  are  physically  realistic:  one 
real  and  one  complex.  Alternatively,  one  may  solve  Eq.  (7)  for  T  and  find  T] 

according  to  Tj  =  r2  +1.  The  proper  k/kQ  is  obtained  from  rj  by  realizing  that 

k  should  have  a  positive  real  part  to  represent  a  wave  propagating  in  the 
postive  x  direction.  Thus  from  the  five  roots  for  T] ,  one  finds  again  five 
roots  for  k/kQ,  two  of  which  have  physical  meaning. 

2.3  Computational  Considerations 


One  may  use  a  complex  (polynomial)  rootfinder  for  the  roots  r  =  iK  of 

2  ^ 

Eq.  (7)  or  Tj  =  (k/kQ)  of  Eq.  (15) ,  and  select  the  physically  acceptable  ones 

by  the  criteria  developed  above.  Alternatively,  one  may  avoid  the 
introduction  of  spurious  roots  that  squaring  caused  by  choosing  the  proper 


branches  in  the  substitution  of  r  according  to  r  =  + 


(k/kQ)  -  1,  again  in 
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accordance  with  the  selection  rules  developed  before.  The  resulting  equations 
are  suitable  to  derive  good  first  approximations  and  iterative  schemes  of 
computation. 


If  one  defines  the  symbol  I  to  indicate  the  positive  real  root  when 

the  argument  is  real  and  positive,  and  the  root  with  the  positive  real  part 
when  the  argument  is  complex  (conform  to  the  sign  convention  in  the  Fortran 

square  root  routine) ,  one  obtains  from  Eq.  (5)  two  equations  for  (k/k^)  :  the 

first  for  complex  k 


i-(k/koy 


1  -  0' 


(k/kD)4] 


+  e/n  =  0 


and  the  second  for  real  k 


(16) 


,  (k/ko)2  -  1  [l  -  n2(k/ko)4]  +  e/n  =  0  .  (17) 

One  necessarily  has  k>kQ  for  real  k,  since  otherwise  the  first  term  of 
Eq.  (17)  would  be  imaginary  and  the  second  term  real,  which  is  impossible. 

The  real  root  for  k/kQ  =  cq/c  may  be  computed  from  Eq.  (17)  by  any  real- 

root  finder.  Simple  iterative  schemes  are  found  as  follows,  starting  from 
previous  knowledge  concerning  the  general  behavior  of  k/k  as  a  function  of 

n  =  u/u  .  ° 

a.  For  0<1 ,  one  finds  that  c  ~  c^  or  Q  (c^/c)  ~  1.  A  good  approximation 

to  Eq.  (17)  is  thus 


(co/c)2  :  Q-‘ 


1  + 


e/n 


In”* 


ij 


1/2  . 


(18) 


and  a  useful  iterative  calculation  to  improve  on  this  first  guess  is 


(CD/C) 


n+1 


=  n 


-1 


1  + 


e/n 


(c  /c)  -  1 

v  o'  'n 


1/2 


,  n  =  1,  2,  3  .  .  .  , 


(19) 


where  one  starts  with  (c^/c) 


2 

1 


given  by  Eq. 


(18). 
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b.  For  fl>l,  one  finds  c  I  cq,  and  thus  a  first  approximation  to  Eq.  (17) 


is 


(c0/c)!  :  1  *  ; 

°  (i-n2)2 

and  an  iterative  scheme  is  given  by 


(20) 


(co/c) 


n+1 


=  1 


'  M*  («o/c)‘J 


n  =  1,  2,  3 


(21) 


'where  the  first  guess  follows  from  Eq.  (20)  . 

Near  0=1  the  convergence  of  these  iterative  schemes  is  very  slow,  and  one 
is  advised  to  use  other  root-finding  schemes;  e.g.,  bisection. 

The  complex  root  of  Eq.  (16)  can  be  found  from  similar  iterative 
computations.  A  good  first  approximation  is  given  by 


(k/k0)!  =  (co/c)!  :  IT1 


1  - 


ie 


nJi-n  * 


1/2 


(22) 


and  by  iteration  one  can  improve  its  accuracy  according  to 


(c  /c)2  1  =  n  1 
v  o'  'n+1 


1  - 


ie 


l-(c  /c) 
v  o'  '  n 


n  =  1,  2,  3 


(23) 


This  scheme  converges  for  all  values  of  fl  but  is  very  slow  close  to  the 
locations  of  the  double  real  roots  discussed  in  the  next  section.  Between 
these  locations  of  double  real  roots  this  latter  iterative  scheme  converges  to 
a  real  root  that  is  physically  not  realistic. 


2.4  Results  for  Fluid-Loaded  Thin  Plate 

The  numerical  examples  in  this  section  refer  to  a  brass  plate,  with  the 
exception  of  Fig.  9.  The  physical  parameters  for  this  material  are  given  in 
Table  1.  Figures  3  through  8  are  given  in  non-dimensional  form  and, 
therefore,  they  depend  only  on  the  loading  parameter  e.  The  thickness  of  the 
plate  enters  only  if  one  wants  to  infer  the  dimensional  frequency  from  the 
value  of  fl.  Also,  in  this  table  derived  quantities  are  given  that  are 
important  to  the  analysis;  e.g.,  the  loading  parameter  e  and  coincidence 
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frequency  w  .  Anticipating  the  results  to  be  given  in  a  later  section  for 

fluid-loaded  constrained  plates,  the  parameters  for  viscoelastic  and 
constraining  layers  used  in  those  examples  are  also  entered  into  Table  1 . 


Figure  3  shows  the  (real  parts  of  the)  roots  of  Eq  (23) ,  expressing 

(k/k^ )  as  a  function  of  the  dimensionless  frequency  ft  =  f/fc  The  validity 

of  the  various  roots  may  be  judged  by  checking  whether  or  not  the  root 
satisfies  Eq,  (17)  when  it  is  real,  and  Eq.  (16)  when  it  is  complex  Valid 
roots  are  designated  by  solid  lines  in  Fig.  3  Notice  that  the  categories 
"realistic”  and  "non-realistic"  do  not  coincide  with  "real”  and  "complex.” 
Curve  1  represents  a  non-radiating  root,  valid  for  the  whole  frequency  range, 
"subsonic,”  since  c/c  =  k  /k<l.  The  value  of  c/c  for  this  root  is  shown  in 

Fig  4,  in  comparison  with  its  value  for  an  unloaded  plate  It  is  interesting 
to  notice  that  the  phase  speed  of  the  flexural  waves  in  the  fluid-loaded  case 
does  not  approach  the  speed  of  bending  waves  in  vacuum,  rather  one  has 


:/c,  -  £  1  ^fi1  ;1  for  fl+0  [5].  For  higher  freque 


asymptotically . 


ncies,  c/cq  approaches  1 


The  solid  part  of  curve  2  in  Fig.  3  represents  the  real  part  of  (k/k^)2 

for  a  radiating  root.  A  close-up  of  curve  2  at  lower  frequencies  is  shown  in 
Fig  5.  (Notice  the  logarithmic  scale  for  the  ordinate  in  Fig.  5.)  One  sees 
that  there  is  a  gap  between  the  two  solid  parts  of  this  curve,  filled  by  a 
pair  of  real  roots  that  are  physically  not  acceptable.  The  branch  points  will 
be  discussed  below. 


Curve  3  of  Fig.  3  shows  the  real  part  of  a  complex  root  pair  that  is 
physically  not  realistic. 

The  value  of  the  non-dimensional  attenuation  a/ k1  corresponding  to  the 
solid  part  of  curve  2,  Fig.  3,  is  shown  in  Fig.  6.  Between  the  two  vertical 
asymptotes  there  is  no  continuation  of  these  complex  roots.  The  situation 
changes  drastically  if  one  introduces  even  a  small  loss  in  the  bending 
stiffness  of  the  plate,  by  setting  D*  =  D(1  +  i^) .  Equation  (16)  now  becomes 


i-(k/kor  [i 


n2(k/k  )4(1  +  117)]  +  e/0  =  0 


(24) 


Now  there  is  one  valid  complex  root  for  all  frequencies,  the  attenuation 
of  which  is  shown  in  Fig.  7.  This  is  the  only  valid  root  left:  the  valid 
real  root  for  the  lossless  case  becomes  complex  and  is  no  longer  physically 
meaningful.  The  curve  in  Fig.  7  shows  a  dip  where  the  relative  attenuation 
approaches  the  value  0.0025.  This  value  follows  from  considering  Eq.  (5)  for 
the  case  of  no  fluid  loading.  Then  its  right-hand  side  equals  zero,  and  one 

has  (k'-ia)4  D(l+ijy)  -  =  0.  Taking  the  imaginary  part  of  this  equation 

one  finds  that  to  first  order  in  the  attenuation  a,  one  has  a/ k1  I  rj/ 4. 

It  is  of  interest  to  study  the  rule  for  the  location  of  the  two  points 
where  the  complex  roots  in  Fig.  5  change  into  a  pair  of  real  roots.  (Remember 
that  the  real  part  of  the  root  depicted  in  the  figure  corresponds  to  a  pair  of 
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complex  conjugate  roots.)  Obviously  at  these  locations  there  is  a  double  real 
root,  found  by  setting  the  derivative  of  the  left-hand  side  of  Eq.  (7)  -with 
respect  to  T  equal  to  0  [10] .  Thus  the  condition  is  that 

5r4  +  6T2  +  1  -  A-2  =  0  ,  (25) 

with  solutions  T2  =  -3/5  _+  (4/25  +  1/5  0  2)1/^2.  Both  double  roots  occur  for 
Q<1,  and  thus  only  the  positive  root  of  the  solution  to  Eq.  (25)  can  be  used. 
The  given  values  for  T  can  be  substituted  into  Eq.  (7) ,  and  the  equation  may 
be  numerically  solved  to  give  two  values  for  fl  and  corresponding  values  for 
c/c^  as  a  function  of  6.  Figures  8a  and  8b  show  the  relevant  curves,  with  a 

few  discrete  points  belonging  to  familiar  materials.  Notice  that  any  material 
with  a  fluid-loading  parameter  €  greater  than  the  right-hand  limit  in  Fig.  8 
will  have  a  continuous  complex  root  for  all  frequencies,  even  without  damping 
in  the  plate.  An  example  is  aluminum  with  e  =  0.39  (see  Fig.  9).  For  a  thin 
plate  loaded  on  both  sides  by  the  same  fluid,  the  right-hand  side  of  Eq.  (2) 
will  be  -2p  instead  of  -p.  This  doubling  of  the  pressure  loading  has  as  a 
consequence  that  the  applicable  e  in  this  case  is  twice  that  of  Eq.  (11) ,  and 
thus  several  of  the  materials  indicated  in  Fig.  8  would  show  this  continuous 
complex  root  for  two-sided  fluid  loading. 

3.  FLEXURAL  WAVES  IN  A  FLUID-LOADED  INFINITE  PLATE,  EXACT  ELASTICITY  THEORY 

The  basic  approach  to  the  description  according  to  exact  elasticity  theory 
of  the  propagation  of  straight-crested  waves  in  a  fluid-loaded  single  plate, 
and  of  reflection  and  transmission  by  such  a  plate,  is  given  elsewhere  [13] . 
The  main  points  are  repeated  here. 

The  dependence  of  the  field  variables  on  the  direction  of  propagation  x 
and  time  t  for  a  harmonic  wave  is  given  by  a  factor  exp [i(wt-kx)] ,  which  is 
subsequently  suppressed.  The  wave  propagation  in  the  plate  is  described  by  a 
potential  for  irrotational  wave  motion 


<f>  =  A  cosh  qz  +  B  sinh  qz  ,  (26) 

s  a 

and  a  potential  for  incompressible  waves 

^  =  C&  cosh  sz  +  Dg  sinh  sz  ,  (27) 

where  q2  =  k2  -k2,  and  s2  =  k2  -  k2  . 

d  s 

Here,  k^  is  the  wavenumber  for  dilatational  waves  with  propagation  speed 

c^,  and  kg  is  the  wavenumber  for  shear  waves  with  speed  cg  in  the  plate 

material.  The  subscripts  a  and  s  in  the  amplitudes  for  the  potentials  (Ag, 

Ba>  and  Dg)  refer  to  antisymmetric  and  symmetric  Lamb  waves,  which  are  the 

generalizations  at  higher  frequencies  for  flexural  and  extensional  waves, 
respectively. 
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One  applies  the  boundary  conditions:  continuity  of  normal  stress  and 
normal  particle  velocity  and  zero  shear  stress  at  the  faces  of  the  plate.  In 
the  case  where  the  same  fluid  is  present  at  both  faces  of  the  plate,  one  has 
in  addition  to  the  pressure  variation  at  the  side  of  the  plate  with  positive  z 
[given  by  Eq.  (3)],  another  pressure  variation  for  negative  z,  given  by 


p_  =  pQ_expi[wt  -  Kxx  +  Kz(z+h/2)]. 


(28) 


The  resulting  matrix  of  the  coefficients  in  the  linear  equations  expressing 
boundary  conditions  is  shown  in  Table  2  The  origin  of  the  coordinate  system 
is  at  the  center  of  the  plate,  and  the  faces  are  at  z  =  +  h/2 

The  various  dispersion  relations  are  obtained  by  equating  the  determinant 
value  of  the  appropriate  matrix  to  zero.  For  fluid  on  both  sides  of  the 
plate,  this  is  the  full  6X6  matrix  of  Table  2;  for  fluid  on  the  positive  side, 
it  is  the  matrix  obtained  by  omitting  row  6  and  column  6  (for  fluid  on  the 
negative  side  one  omits  row  5  and  column  5;  but,  of  course,  the  result  is  the 
same) .  For  a  plate  in  vacuum,  rows  5  and  6  and  columns  5  and  6  are 
suppressed. 


To  appreciate  the  effects  of  using  exact  elasticity  theory,  the  phase 
speeds  of  flexural  waves  and  extensional  waves  are  shown  in  Fig  10  and  are 
compared  with  the  results  of  thin-plate  theory .  The  results  are  given  for  a 
10-cm  thick  brass  plate;  the  parameters  are  listed  in  Table  1. 


One  sees  that  the  value  of  the  frequency  for  which  the  flexural  wave  speed 
is  equal  to  the  sound  speed  in  the  medium,  water,  is  shifted  from  3218  Hz  for 
thin-plate  theory  to  5199  Hz  for  exact  elasticity  theory. 


The  phase  speed  for  a  plate  loaded  on  one  side  is  shown  in  Fig.  11 
Notice  that  in  exact  elasticity  theory  there  are  only  two  roots  in  the 


k2  - 

o 


In  thin- 


dispersion  equation,  depending  on  the  sign  of  K  =  + 

z  — 

plate  theory  there  are  five  roots  of  the  polynomial  Eq.  (7) .  The  appearance 

of  this  figure  shows  features  similar  to  those  encountered  in  the  case  of 

thin-plate  theory,  see  Figs.  4  and  5.  There  is  again  a  non-radiating  root  at 

all  frequencies  with  asymptotic  value  c  for  i*«:  the  radiative  root  consists 

o 

of  two  branches ,  interrupted  by  a  frequency  range  where  there  are  two  real 
roots  that  are  not  physically  meaningful.  The  points  where  this  transition 
takes  place  are  at  frequencies  875  and  3635  Hz.  The  corresponding  curve  for 
a/k'  is  shown  in  Fig.  12. 


4.  HYBRID  MODEL  OF  FLUID-LOADED  CONSTRAINED  PLATE 


4.1  Theory 

A  model  has  been  developed  [3]  for  constrained  layer  damping,  whereby  the 
base  plate  is  described  by  exact  elasticity  theory,  while  the  constraining  and 
damping  layers  are  treated  according  to  the  Kerwin  model  [1] . 
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This  hybrid  model  may  be  easily  extended  to  include  fluid  loading  by 
combining  the  modified  boundary  condition  due  to  the  two  extra  layers  with  the 
fluid  loading  implicit  in  the  matrix  of  Table  2. 

As  derived  in  Ref.  13.  the  shear  stress  r  at  the  side  of  the  plate 

adjacent  to  the  elastomer  layer  is  not  zero,  as  for  an  unconstrained  plate, 
but  is  given  by 


1  -  (c/c3)2  G2u 

TzX  "  1  -  (c/c3)2+g  h2  ' 

where  c^  is  the  thin-plate  extensional  wave  speed  of  the  constraining  layer. 

The  elastomer  layer  has  a  complex  shear  modulus  and  thickness  h2 .  The 

tangential  displacement  of  the  surface  of  the  plate  is  u.  The  shear  parameter 
g  is  given  by 


(29) 


g  =  G2(l-^)/(k2  E3h2h3) 


(30) 


where  E^  is  Young’s  modulus,  Poisson’s  ratio,  and  h^  the  thickness  of  the 
constraining  layer. 


The  displacement  u  in  Eq.  (29)  may  be  expressed  in  terms  of  the  amplitudes 

A  ,  B  ,  C  ,  D  of  Eqs.  (26)  and  (27).  The  result  of  inserting  these  into  the 
s  a  a  s 

expression  for  the  shear  stress  Eq.  (29)  is  a  modification  of  the  first  four 
elements  of  the  second  row  of  the  matrix  in  Table  2,  assuming  that  the 
constraining  layer  is  at  the  positive  z  side  of  the  plate.  The  modified 
matrix  elements  m2 ^  are  (using  the  half  thickness  d^  =  h  /2  of  the  base  plate) 

2 

m2J  =  2ikqdi  cosh(qdl  )  +  iFkdj sinh(qd1 )  , 


m22  =  (k2 +s2  )d2  cosh(sdl  )  +  Fsdj  sinh(sdl  )  , 

(31) 

m  =  2ikqd2 sinh(qd  )  +  iFkd  cosh (qd,  )  , 

l  o  1  1  1  1 

m24  =  (k2 +s2  )d2  sinh(sd1 )  +  Fsd^  cosh(sdl )  , 


■where  F  -  0.5  (G2  /Gj )  f\ )  [1  “  (c/c3)2]/[l  ~  (c/c3)2  +  g]  .  G1  is  the  shear 
modulus  of  the  base  plate  material. 

The  dispersion  relation  for  the  case  where  constraining  layer  and  fluid 
are  on  the  same  side  is  obtained  by  suppressing  row  6  and  column  6  in  the 
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matrix  If  constraining  layer  and  fluid  loading  are  on  opposite  sides,  column 
5  and  row  5  are  suppressed  If  both  sides  are  loaded  by  fluid,  the  full  6X6 
matrix  is  used 

4.2  Evaluation  of  Damping  Efficiency  in  Fluid-Loaded  Constrained  Plate. 

The  purpose  of  constrained-layer  damping  is  to  absorb  acoustic  energy 
present  in  a  vibrating  plate  or  bar  by  means  of  a  viscoelastic  layer. 
Therefore,  one  tries  to  maximize  the  energy  absorbed  by  studying  the 
attenuation  coefficient  for  viscous  damping  a as  a  function  of  the  parameters 

of  the  experimental  condition.  When  a  constrained  layer  is  loaded  by  fluid, 
radiation  will  take  place  at  certain  frequency  ranges  and  it  is  obviously 
important  to  assess  what  fraction  of  the  total  attenuation  is  due  to  loss  by 
radiation.  Apart  from  the  attenuation  by  viscous  damping,  one  should 
therefore  also  consider  the  damping  efficiency,  defined  as  the  fraction  of  the 
total  absorbed  energy  that  is  lost  to  viscoelastic  dissipation. 

Let  W  be  the  power  propagated  in  the  straight-crested  wave  per  unit  of 
width  of  the  plate.  Then  this  damping  efficiency  may  be  computed  in  the  case 
of  the  hybrid  model  by  evaluating  the  rate  of  energy  loss  by  radiation, 
(dW/dx)r,  and  the  rate  of  energy  loss  due  to  viscoelastic  absorption, 

(dW/dx)y,  per  unit  of  distance  along  the  plate. 

The  energy  loss  by  radiation  follows  from  evaluating  the  intensity  vector 
1,  expressed  in  terms  of  the  amplitudes  of  particle  velocity  vector  v  and 
pressure  p  by  I  =  -Re(p  v) ,  where  Re  indicates  the  real  part  of  a  complex 
quantity,  and  p  is  the  complex  conjugate  of  p.  The  pressure  at  the  positive 
side  of  the  plate  p  is  given  by 


P+  =  PQ+expi[wt  -  kx  -  Kz(z-h/2)], 

where  k  =  k1  -  ia,  and  K  =  K1  +  i  K" .  The  velocity  vector  v  is 

z  z  z 


(32) 


V 

+ 


^oCo 


k1  -ia 

k 

o 


K'+iK" 
z  z 

k 

o 


The  intensity  vector  is  then  given  by 


(33) 


!  = 


I  |2 

lP+ol 

,  ,  K' 
k'  z 

2  p  c 
ro  o 

k  ’  k 
o  oJ 

and  its  absolute  value  is  given  by 


(34) 
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I  = 


CN 

2 

1/2 

Kl2  . 

+ 

K' 

zJ 

2  p  c 
ro  o 


i  i2 

2P0cf 


(35) 


inhere  is  the  phase  speed  of  the  radiated  wave  introduced  before.  The 

expression  Eq.  (35)  has  the  well-known  form  for  plane  waves  of  the  free  wave 
speed  cq  is  replaced  by  the  speed  c^  for  waves  influenced  by  the  presence  of 

the  plate.  The  energy-loss  rate  per  unit  of  distance  due  to  radiation  is  thus 

(suppressing  a  factor  e  2ax) 


dW 


dx 


=  I  cos 8  =  — 


i  ipr 


cos5  =  t: 


P 

2 

K' 

1  ro+ 

z 

o 

o 

D 

k 

L  O  > 

2  p  c. 
ro  f 

The  computation  for  the  other  side  gives  a  similar  result, 


(36) 


dw 

i 

Po- 

dx 

"  2 

T* 

Pnc, 

cos<5 


(37) 


The  rate  of  energy  loss  per  unit  of  distance  in  the  x  direction  in  the 
viscoelastic  layer,  according  to  Kerwin’s  model,  is  given  by 


dW 

dx 


v 


(38) 


where  £  is  the  amplitude  of  the  shear  displacement  in  the  elastomer  layer, 
equal  to  the  difference  in  the  displacement  of  the  top  and  bottom  of  the 
layer.  This  may  be  written  as 


l£|2  =  ie/ul2  lu/c^l2  ,  (39) 

where 

1  -  (c/c )! 

£/”  - - h -  •  (40) 

1  -  (c/c3)2  +  g 

The  field  quantities  PQ+>  PQ_,  anc^  (C/u)  can  5e  obtained  from  the  matrix  in 
Table  1.  One  can  then  find  the  damping  efficiency  from 
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! dW/dx I 

D  =  I  dW/dx  I  r+ !  dW/dxl 


(41) 


4.3  Comparison  of  Damping  Efficiency  in  Various  Models 

The  concept  of  damping  efficiency  was  defined  in  section  4.2  as  the 
fraction  of  the  total  energy  loss  that  is  due  to  damping  by  the  viscoelastic 
layer.  It  is  quite  common  to  evaluate  this  efficiency  by  computing  it  from 
the  separate  cases  of  a  single  fluid-loaded  plate  and  a  constrained  plate  in 
vacuum  This  implies  the  tacit  assumption  that  the  radiation  by  a  fluid- 
loaded  constrained  plate  is  the  same  as  the  radiation  by  a  simple  fluid-loaded 
plate,  and  that  the  viscoelastic  damping  in  a  fluid-loaded  constrained  plate 
is  the  same  as  that  in  a  constrained  plate  in  vacuum  It  ignores  a  possible 
interaction  between  the  fluid-loading  and  constraining  layers  in  the  composite 
situation. 

In  this  section,  the  validity  of  this  assumption  is  investigated  by 
comparing  the  damping  efficiency  for  a  fluid-loaded  constrained  plate  with  a 
linear  combination  of  the  separate  cases  of  a  single  plate  loaded  by  fluid  and 
a  constrained  plate  in  vacuum  The  damping  efficiency  for  the  first  case  is 
computed  according  to  the  expression  Eq.  (41) .  For  the  linear  combination  the 
attenuation  constant  for  the  single  fluid-loaded  plate  a  is  combined  with  the 
attenuation  constant  a  for  the  constrained  layer  in  vacuum  to  give  the  total 
attenuation  by 


aT  ~  “v  +  ar  ’ 


(42) 


and  damping  efficiency  by 


€r,  = 


a  +a 
v  r 


(43) 


In  the  following  numerical  examples  the  base  plate  is  10-cm-thick  brass.  The 
elastomer  layer  is  1.24-mm  thick,  its  shear  modulus  is  chosen  to  be  constant 
for  the  whole  frequency  range.  Of  course,  this  is  physically  unrealistic 
since  the  property  of  viscoelasticity  is  by  nature  dependent  on  frequency. 

The  constancy  of  G,  though,  brings  out  the  geometry  and  frequency  dependence 
of  the  damping,  separated  from  viscoelastic  behavior  as  a  function  of 
frequency.  A  value  of  G2  =  10(l+i)  MPa  was  chosen  The  constraining  layer  is 

2.48-mm  thick.  Table  1  shows  the  parameters  of  the  constrained  plate  and  of 
the  acoustic  medium. 

Figure  13  shows  the  comparison  for  a  constrained  plate  loaded  by  a  fluid 
at  the  side  of  the  constraining  layer  with  the  separate  cases  of  a  single 
plate  with  a  fluid  layer  and  a  constrained  layer  in  vacuum.  One  sees  clearly 
the  two  ranges  of  high  damping  due  to  radiation  and,  in  between  these,  the 
hump  due  to  viscoelastic  damping.  It  appears,  at  least  at  the  scale  of  the 
graph,  that  the  total  attenuation  is  practically  due  to  one  or  the  other 
effect  with  only  some  visible  discrepancy  at  the  higher  end  of  the 
viscoelastic  hump. 
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The  comparison  of  the  damping  efficiency  computed  for  the  two  cases  (the 
first  being  the  linear  combination  of  a  single  plate  in  the  fluid  with  a 
constrained  plate  in  vacuum,  the  second  being  a  fluid-loaded  constrained 
plate)  is  shown  in  Fig.  14.  Again,  for  a  rather  large  frequency  range  the 
differences  are  small.  At  the  lower  frequency  end  there  is  a  discrepancy  that 
apparently  corresponds  to  the  fact  that  the  absorption  in  the  viscoelastic 
layer  is  influenced  by  the  presence  of  the  fluid.  At  the  high-frequency  side 
one  sees  that  the  efficiency  computed  from  the  linear  combination  turns  up 
sharply  at  about  40  kHz.  This  is  due  to  an  effect  already  discovered  [13] . 

For  high  frequencies  the  loading  of  one  face  of  the  plate  by  a  fluid  forces 
the  action  (i.e.,  the  displacement  of  the  plate  particles)  to  the  unloaded 
side  so  that  the  radiation  drastically  diminishes.  This  suggests  that  in  this 
frequency  range  it  is  better  to  damp  a  wave  by  applying  the  constraining 
layers  at  the  side  of  the  plate  opposite  to  the  loading  fluid.  The  effect  of 
this  can  be  seen  in  the  total  attenuation  (Fig.  15)  where  attenuation  due  to 
radiation  falls  off  for  both  cases  (fluid  on  the  same  side  or  opposite  side  to 
the  constraining  layers)  while  the  viscoelastic  damping  remains  intact  for  the 
latter  case.  The  effect  on  damping  efficiency  is  clearly  visible  in  Fig.  16. 

When  the  attenuation  for  one-sided  loading  is  compared  with  the  case  of 
fluid-loading  on  both  sides  (Fig.  17),  one  sees  that  indeed  the  radiation 
stays  high  for  the  symmetric  situation.  Accordingly,  the  damping  efficiency 
(Fig.  18)  shows  a  steady  decrease  with  f  and  is  well  approximated  by  a  linear 
combination  of  separate  radiation  and  viscoelastic  damping  (dashed  curve) . 
Notice  that  in  Figs.  16  through  18  only  the  high-frequency  part  of  the  curves 
is  shown.  The  scale  of  the  abscissa  is  kept  the  same  to  facilitate  comparison 
with  Figs.  12  through  14. 


5.  CONCLUSIONS 

This  study  shows  that  for  a  sizable  frequency-range,  it  is  a  good 
approximation  to  obtain  the  damping  and  radiation  in  a  fluid-loaded  composite 
plate  by  a  linear  combination  of  the  attenuation  in  a  single  plate  in  a  fluid 
and  a  composite  plate  in  vacuum. 

At  both  the  high-  and  low-frequency  ends  of  this  range,  though, 
discrepancies  occur  between  the  linear  combination  of  the  two  cases  of  a 
single  plate  in  a  fluid  and  a  constrained  plate  in  vacuum  and  the  computation 
based  on  a  fluid-loaded  constrained  plate  by  the  hybrid  model.  This 
necessitates,  in  general,  the  study  of  the  complete  fluid-loaded  composite 
situation,  in  order  to  predict  the  total  attenuation,  and  especially  the 
damping  efficiency. 
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Table  1.  Material  Parameters  and  Derived  Quantities 


1.  Base  plate,  brass 
thickness 
Young’s  modulus 
Poisson’s  ratio 
density 


h,  = 


E,  = 


z/,  = 


10  cm 
104  GPa 
0.37 


p  =  8500  kg/m 


shear  modulus 

extensional  wave  speed 
(thin-plate) 


shear  wave  speed 
bending  stiffness 
fluid-loading  parameter 
coincidence  frequency 
(thin-plate) 


c 

D5 

e 

f 


2.  Viscoelastic  layer  (hypothetical) 
thickness 
shear  modulus 

loss  tangent  in  ()  =  1.0 
shear  modulus 


hg  =  1.24  mm 
G2  =  10  MPa 


3.  Constraining  layer,  aluminum 
thickness 
Young’s  modulus 
Poisson’s  ratio 
density 


hg  =  2.48  mm 
E3  =  71  GPa 
u3  =  0.33 


4.  Fluid,  fresh  water 
density 
wave  speed 


p 3  =  2700  kg/m 


PQ  =  998  kg/m3 
c  =  1481  m/s 

r\  t 


=  38.0  GPa 
=  3765  m/s 


=  2113  m/s 
'=  10.0  MN  m 
=  0.086 
=  3212  Hz 
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Table  2  -  Matrix  of  Coefficients  of  Equations  Describing  Wave  Propagation  in  Fluid-Loaded  Single  Plate. 

Exact  elasticity  theory.  Here  d  is  the  half  thickness  of  the  plate  d  ~  h/2 
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qd  cosh(qd)  -ikd  cosh(sd)  -qd  sinli(qd)  ikd  sinh(Bd) 
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CONSTRAINING 

LAYER 


VISCOELASTIC 

LAYER 


BASE  PLATE 


Fig*  1  -  Geometry  of  conatrained  layer* 


Fig.  2 


-  Relation  between  ware- amplitude 


WAVE  IN  FLUID 

X 

WAVE  IN  PLATE 

variation  in  the  plate  and  in  the  fluid. 
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Fig.  8  -  Beal  part  of  the  roots  TJ  =  (k/k  of  Eq. 

o 

plate.  The  solid  (parts  of)  curves  are  physically 


Curve  1  Real  subsonic  root 
Curve  2  "  radiating  root 

- . —  two  real  roots 

The  branch  point  is  located  at  f/f 
Curve  8  - complex  root 


22  for  flexural  waves  in  fluid-loaded 
realistic. 


=  0.726,  (k/k  )2 
o 


1.136 


Fig.  4  -  Relative  ware  speed  of  non-dissipative  subsonic  flexural  ware  in  thin  plate  as 
function  of  relative  frequency.  Loaded  by  fluid  on  one  side.  _  in  Ta 
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Fig.  6  Fluid  loaded  thin  plate.  Acceptable  complex  root.  -  Non-realistic 

real  root*.  Extension  of  curre  2,  Fig.  3  for  low  frequency.  The  branch  points  are  located 

2  o 

*/*  =  0*725}  (k/k  )  =  1.136  and  f/f  =  0.0276;  (k/k  )  =  16.69. 

Co  CO 


Fig.  6  -  Fluid-loaded  thin  plate.  Non-dimensional  attenuation  a/ k1 
for  radiating  root  of  Fig.  6 
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Fig,  7  -  Fluid-loaded  thin  plate.  Non-dimensional  attenuation,  OL/ k1  for  plate  with 
loss  tangent  Tj  =  0.01  in  bending  stiffness  D. 


Fig.  8  -  Thin  plate  loaded  on  one  side  by  water, 
of  fluid-loading  parameter  €. 

Discrete  points t  1  -  lead  4  -  nickel 

2  -  silver  6  -  steel 

3  -  brass 


Values  of  double  real  roots,  as  a  function 

a)  ft  =  U/UI 

c 

b)  c/c 

o 
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Fig.  9  Thin  aluminum  plate  loaded  on  one  side  b7  water.  Non-dimensional  attenuation 
a/ k'  as  &  function  of  relatire  frequency  f/f  . 


T-00  too 

FREQUENCY  (kHz) 


KQQ 


10 

? 

1 

£  S.0 

ft 

l 

1.0 


ol_ 
o  so 


Fig.  10  -  Phase  speed  of  waxes  in  infinite  brass  plate  of  10  cm  thickness  in  racuum 
Antisymmetric  (flexural)  ware 
1 •  Exact  theory 

2.  Thin-plate  theory 
Symmetric  (exteneional)  ware 

3.  Exact  theory 

4.  Thin-plate  theory 
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Fig.  IX  Phase  speed  of  10-cm  thick  brass  plate,  loaded  by  water  on  one  side.  The  solid 
curves  represent  physically  acceptable  roots. 

X.  (solid  parts)  radiative  root 
2.  Non-dissipative  root 


Pig.  12  -  Non-dimensional  attenuation  Q./ k1  of  radiating  wares  in  10-cm  thick  brass  plate 
loaded  on  one  side  by  water. 


24 


NRL  MEMORANDUM  REPORT  6067 


Fig.  18  -  Comparison  of  attenuation  coefficient*.  -  Total  dimensionless 

attenuation  a/ k'  of  constrained  plate,  loaded  by  water  on  the  side  of  the  damping 

layers.  -  Attenuation  of  single  plate  loaded  by  water  on  one  side. 

“  *-*-*-*  Damping  of  constrained  plate  in  Ttcuum - 


Fig.  14  -  Comparison  of  damping  efficiency.  -  Constrained  brass  plate  loaded  by  water 

on  side  of  damping  layers.  -  Computed  by  linear  combination  of  radiation  damping 

of  single  plate  loaded  on  one  side,  and  viscoelastic  damping  of  constrained  plate  in  vacuum. 
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Fig.  16  -  Comparison  of  total  attenuation.  - - -  Constrained  bras,  plate  loaded  by  water 

on  side  of  damping  layers.  -  Constrained  bras,  plate  loaded  by  water  on  opposite 

side  from  damping  layers* 


Fig.  16  -  Comparison  of  damping  efficiency.  -  Constrained  layer  loaded" by  water  on 

side  of  damping  layer.  -  Constrained  layer  loaded  by  water  on  opposite  side 

from  damping  layers • 
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Fig.  17  -  Comparison  of  total  attenuation.  -  Constrained  brass  plate  loaded  by  water 

on  side  of  damping  layers.  -  Constrained  brass  plate  loaded  by  water  on  both  sides. 


Fig.  18  -  Comparison  of  damping  efficiency.  -  Constrained  layer  loaded'by  water  on 

both  sides.  -  Computed  by  linear  combination  of  attenuation  of  radiation  damping  for 

single  plate,  loaded  on  both  sides,  and  riscoelastic  damping  of  constrained  plate  in  vacuum. 


27 


U 230623  I 


3  o 

<  -n 
■n  2 

§2 

*  > 

I: 

<  2 

5  w 


c  m 
w  g 

rn  C/5 

«/» 

CO 

O 

O 


$ 

02 

z 

> 

< 

V» 

> 

3* 

f— 

5* 

CD 

X 

m 

O 

C/5 

3 

m 

> 

o 

X 

p 

o 

X 

N3 

r~ 

O 

> 

CJ 

U3 

•vJ 

Ul 

O 

X 

ui 

> 

O 

H 

O 

O 

o 

X 

-< 

O 

m 

-o 

> 

x 

H 

2 

m 

Z 

H 

O 

m 

H 

X 

m 

Z 

> 

< 

■< 


o 

m  -o 

3  ° 

>  05 

X  H 

-i  > 

2  £ 
n  m  rn 

p  ?  Z  > 

C  o  H  z 

w  w  O  ° 

C/5  -*  ^  -n 
_  o>  h  m 

2  x  K 

>  m  c/> 

F  2  2 

>  £ 

<  o 


X 

X 

o 


